Commensurability oscillations due to pinned and drifting orbits in a two-dimensional 

lateral surface superlattice 
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We have simulated conduction in a two-dimensional elec- 
tron gas subject to a weak two-dimensional periodic potential, 
VxCos{2-Kx/a) + Vy cos{2ny / a) . The usual commensurability 
oscillations in pxx{B) are seen with Vx alone. An increase of 
Vy suppresses these oscillations, rather than introducing the 
additional oscillations in pyy{B) expected from previous per- 
turbation theories. We show that this behavior arises from 
drift of the guiding center of cyclotron motion along contours 
of an effective potential. Periodic modulation in the magnetic 
field can be treated in the same way. 



The behavior of electrons in a periodic potential lies 
at the heart of solid state physics and continues to yield 
surprises. Motion in a controllable one (ID) or two- 
dimensional (2D) potential can be studied with a lateral 
surface superlattice (LSSL). The electrons typically lie in 
a high-mobility 2D gas in a semiconducting heterostruc- 
ture, and the periodic potential is applied through an 
array of metal gates whose bias can be varied. Alterna- 
tively a patterned stressor may be used, in which case 
the dominant potential is piezoelectric; this has a lower 
symmetry than the stressor, which will prove important. 

The aim of using LSSLs is often to explore quantum- 
mechanical effects, such as Bloch oscillation and the Hof- 
stadter butterfly, but the period of the potential is too 
long in most current devices. Instead, the dominant ef- 
fects seen in ID LSSLs are commensurability oscillations 
(COs) in the magnetoresistance B. These can be ex- 
plained semiclassically B from interference between cy- 
clotron motion and the periodic potential. Consider a 
sinusoidal potential energy, V{x) — T4 cos(27rx/a). The 
interference causes a drift along the equipotentials, par- 
allel to the y axis, which contributes to the conductivity 
ayy and the resistivity pxx- 
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Here Jq is a Bessel function of the first kind, po is the 
resistivity at i? = 0, Hs the mean free path, i?^ = vp/ujc 
is the cyclotron radius, lOc ~ eB/m is the cyclotron fre- 
quency, E-p is the Fermi energy and v-p the Fermi velocity. 
No effect on pyy is expected in this approach. Quantum- 



mechanical analysis I3M| yields a similar result but with 
small contributions to pyy. Overall agreement between 
theory and experiments on ID LSSLs is excellent, even 
for the strong piezoelectric potentials in strained LSSLs 



Now consider a simple 2D potential energy, 
V{x, y) = Vx cos(27rx/a) -f Vy cos{2TTy/a). 
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To avoid ambiguity we take Vx > Vy. An extension of the 
semiclassical theory [Q shows that Vx continues to gener- 
ate oscillations in pxx{B) according to Eq. (||), and Vy has 
the same effect on pyy{B). In contrast to ID potentials, 
there is little confirmation of this plausible result. An 
early experiment g used a holographic technique 0] in 
two steps. A ID grating was produced first, and showed 
strong COs in the longitudinal resistivity as expected. 
The sample was next illuminated with an orthogonal ID 
pattern to produce a 2D grid. However, this combined 
pattern did not produce similar, strong COs in both pxx 
and Pyy , as expected from the extension to the semiclas- 
sical model; much weaker oscillations with the opposite 
phase were seen instead. Many subsequent measurements 
have used different modulation techniques. Virtually the 
only common feature between them is that the COs are 
generally weak, confirming this most significant feature 
of the holographic experiment. 

We have performed simulations of conduction in a 2D 
LSSL to address this issue, and find that Vx and Vy do 
not contribute independently. Instead, the introduction 
of Vy suppresses the oscillations in pxx rather than in- 
ducing oscillations in pyy. We explain this with a simple 
picture based on drift of the guiding center of cyclotron 
motion along contours of an effective potential. Trajec- 
tories can drift or be pinned, and the latter suppress 
the magnetoresistance. Remarkable behavior is found 
when higher Fourier components dominate this poten- 
tial, which should be detectable in experiments. 
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FIG. 1. Simulation of diagonal components of resistivity 
tensor for square two-dimensional superlattices with period 
200 nm, T4 = ImeV, and several values of Vy. Calculated 
points are joined by lines for clarity and the thick curve shows 
the semiclassical result (Eq. Ill) for Vy — 0. Two simulations 
are plotted for Vy = Vx and lie on top of each other. The 
inset curve shows the estimate (Eq. 0) of the effect of Vy on the 
oscillations due to V^ alone, with points from the simulations. 

To simulate conduction we solved the classical equa- 
tions of motion for electrons moving in the potential en- 
ergy given by Eq. (g) and a normal magnetic field, with 
a constant probability of isotropic scattering. The super- 
lattice had period a = 200 nm in GaAs with 3 x 10^^ m^^ 
electrons of mobility 50m^ V~^ s~^. The resistivity ten- 
sor was deduced from the velocity autocorrelation func- 
tion and its diagonal elements are plotted as a function 
of the magnetic field in Fig. |^. We held V^ = 1 meV 



and raised Vy from zero to V^ 



In the ID limit, Vy = 0, 



the usual oscillations are seen in pxx with no structure in 
Pyy^ in excellent agreement with Eq. (|l|). Recall that the 
existing theory predicts that an increase of Vy from zero 
should induce oscillations in pyy without affecting pxx- 
Instead, we see no oscillations in pyy while those in pxx 
are suppressed. Most strikingly, there are no commensu- 
rability oscillations at all in the symmetric 2D limit where 
Vy = Vx- The only structure that remains is a positive 
magnetoresistance at low fields. This arises from mag- 
netic breakdown |^,|lO| and it has recently been suggested 
PI that successive breakdowns may lead to further os- 
cillations of quantum-mechanical origin. 




FIG. 2. (a) Drifting and (b) pinned trajectories in a 
square superlattice with period 200 nm shown by the grid, 
14 = ImeV, Vy = |meV, and B = 0.72 T. (c) Chaotic 
trajectory for B = 0.28 T; note the different scale. 

An explanation of this behavior follows from the tra- 
jectories taken by electrons in the simulation. Typical ex- 
amples starting from different points are shown in Figs. 
|(a) and (b) for 14 = ImeV and Vy = i meV. The 
magnetic field B = 0.72 T, corresponding to the largest 
peak in Fig. |l|. There is no scattering and the trajecto- 
ries run for 100 ps, considerably longer than the lifetime 
T = 19 ps if scattering had been included. There is no 
sign of the chaos seen in weaker magnetic fields |12] . The 
underlying motion is clearly a cyclotron orbit and the 
overall trajectories can be divided into two classes. Fig. 
H(a) shows the cyclotron orbit drifting along y. This is 
perpendicular to the wavevector of the stronger poten- 
tial component, and is the only type of trajectory seen 
in the ID limit, Vy ^ 0. Such motion contributes to 
ayy and pxx- No electrons were found to drift along x 
and we therefore expect no effect on axx and pyy. Tra- 
jectories of the second class are pinned, as in Fig. @(b). 
The cyclotron orbit is distorted and precesses but makes 
no overall displacement. Such orbits make no contribu- 
tion to conduction in the limit of large u>cT and therefore 
suppress the magnetoresistance. All trajectories become 
pinned in the symmetric 2D limit, Vy = Vx, quenching 
the COs. This analysis of the trajectories therefore shows 
that COs are reduced in magnitude because of pinned or- 
bits, and are seen only in pxx if Vx > Vy 
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FIG. 3. Contour plots of effective potential energy for 



a square superlattice. (a) Vy 



K" 



with all contours 



closed, (b) V;; 



2 ^ 



showing two bands of open con- 



tours (shaded) running parallel to the Y axis. 

The trajectories are only weakly distorted from regular 
cyclotron motion. We therefore focus on motion of the 
guiding center B , which drifts at a velocity given by 

v(^) {X, Y) = S/V'^iX, Y) X B/eB^. (3) 

The effective potential energy V^ {X, Y) is the periodic 
potential energy (Eq. y) averaged over the perimeter of 
a cyclotron orbit centered on {X,Y), which reduces T4 
and Vy equally by a factor of Jo(27ri?c/a). This depends 
on magnetic field through the cyclotron radius Re- Eq. 
(0) shows that the guiding center drifts along contours 
of V''^{X,Y). Two examples are plotted in Fig. |. All 
contours are closed for a symmetric effective potential 
with V^^ — Vy^ [Fig. 3(a)]. All trajectories are there- 
fore pinned as in Fig. 2(b). Fig. p|b) shows the effect 
of breaking the symmetry with Vy = ^V^^ . This in- 
troduces a band of open contours, shaded in the plot, 
running parallel to the Y axis. The guiding center can 
drift along these and the deviations of the contours from 
straight lines leads to the lateral oscillations seen in Fig. 
|(a). 

The change in conductivity can be estimated [p]p3| 
from 
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It is assumed that only drifting orbits need be considered. 
We first find the average drift velocity v^'^) for each orbit, 



(d)l2 :, 



which leaves only Vy in our case. The square [vy 
then averaged over all orbits in the unit cell to give ^(Jyy . 
A difficulty is that Eq. (0) is valid only if the lifetime r is 
much larger than the periods of the drifting orbits. This 
always fails for trajectories on the boundary of the open 
region because they go through stagnation points in the 
middle of each edge of the unit cell, but is easier to satisfy 
for the majority of orbits. 

The open orbits are complicated and we therefore make 
several approximations to estimate Eq. (Q). Start from 
the ID limit, Vy = 0, in which case all orbits drift with 

v^y'^\x) = {2TTV^^/eBa)sm{2TTX/a). The introduction 
of Vy affects this in two ways. First, the fraction of the 
unit cell occupied by drifting orbits is reduced to 
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We replace the areas of drifting orbits shown in Fig. ^(b) 
by bands along Y of the same area centered on X = jU 

and |a, and average [vy {X)]'^ over the remaining area. 
The second effect of Vy is to make the drifting orbits 
sinuous, which reduces their average velocity along Y. 
The most rapid orbit is through the symmetry point 
{j,j)a. Its period is increased by a factor of {2/'k)K 
compared with Vy — 0, where K is the complete elliptic 
integral of first kind jlj] with modulus k — VyjVx- We 
apply this factor to all orbits. Combining the two effects 
lead to the approximate resistivity 
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This is plotted in the inset to Fig. |^. It reduces correctly 
to the symmetric and ID limits, and agrees well with the 
simulations. 

These results show that COs are much harder to ob- 
serve in 2D potentials because the symmetry must be 
broken. They are also less robust in 2D. If V^ = the 
guiding center drifts along Y , which does not change the 
potential seen by the electron. In 2D, however, the poten- 
tial changes and the picture based on the guiding center 
will be valid only if its drift during one cyclotron period 
is much smaller than the unit cell. This leads to the 
condition eBa '3> 2'K{mV°^Y/'^ for COs to exist, which 
is similar to that for normal diffusion rather than chaos 
in a much stronger, symmetric potential [n2|. Using the 
envelope of the Bessel function to relate V^ and Vx al- 
lows this to be rewritten as VF{eBa)^ » m(27rT4)^. This 
becomes B :^ 0.3 T for the conditions used in Fig. ||, 
which is in reasonable agreement with the simulations. 
We would also expect the COs to become more robust 
as Vy is reduced, which is seen. Motion becomes chaotic 
when the condition on B is violated and a typical trajec- 
tory is shown in Fig. |(c) for B = 0.28 T. 

A general feature that follows from the drift of the 
guiding center along contours, Eq. (||), is that there can 
be only one average direction of drift. By symmetry this 
must be along AT or F if only Vx and Vy are present, 
so it is impossible to have oscillations in both pxx and 
Pyy. New features appear when further Fourier com- 
ponents are added to the potential. Consider the sim- 
plest 'diagonal' component, Vi,i cos[27r(a; -I- 2/)/a]. The 
average over a cyclotron orbit gives an effective poten- 
tial energy Fff = l^i,i Jo(2y27ri?c/a). This has a differ- 
ent dependence on magnetic field from Vx and Vy, and 
Vl^i can therefore dominate the behavior near zeroes of 
J{){2iTRc/a) where the fundamental terms vanish. 




FIG. 4. Effect of diagonal component Vj^i on contour plots 
of effective potential energy for a square superlattice with 
Vy = ^V^ . (a), (b) Region of open contours distorts as 
Vi'^i is raised but drift remains parallel to Y on average. 

(c) Open contours vanish for Vi^ = V^^; all orbits are pinned. 

(d) Open contours reappear for Vii > V^ but now run di- 
agonally. 

The effect of raising Vi^ from zero is displayed in Fig. 
^, holding Vy^ — ^V^^ . The region of open contours 
parallel to Y is distorted and shrinks as V^f rises. It col- 
lapses when y/^f = V^^, all orbits are pinned, and COs 
are quenched. Open contours reappear when V]°f > V°^ 
but are now parallel to the diagonal Y = —X. This di- 
agonal drift induces a peak in both pxx and pyy instead 
of the expected minimum in p^x- Such a mechanism may 
contribute to the antiphase effects seen in some experi- 
ments [pi , particularly if the fundamental potential com- 
ponents are well balanced, and the corresponding COs 
are suppressed. Note, however, that no COs of any pe- 
riod will be seen within this model unless there is some 
asymmetry in the potentials. 

Commensurability oscillations can also be induced by 
a 2D periodic magnetic field, Bz — Bq + Bm{x,y) 
|l3| , |l5|Jl6| . This can be analyzed in the same way with 
an effective potential energy given by V^{X,Y) — 
EfB^{X,Y)/Bo. Here the magnetic field B^{X,Y) 
is averaged over the area of a cyclotron orbit centered 
on {X,Y), rather than its perimeter. This introduces a 
factor of 2Ji{9)/9 with 9 = 2-KRc/a for the fundamen- 
tal Fourier components, rather than Jo{9). Experiments 
show that COs in a square array of magnetic elements 
appear only in the direction of in-plane magnetization 
[pTI in accord with contours of B^ presented there. 

We have shown that the commensurability oscillations 
in a two-dimensional superlattice are quite different from 
the superposition of one-dimensional results. The addi- 



tion of the potential energy Vy cos{2iTy/a) suppresses the 
oscillations in Pxx{B) due to Vx coa{2iTx/a) alone, rather 
than adding new oscillations in pyy{B). There are no 
oscillations at all in a symmetric potential, Vx = Vy, pro- 
vided that only open orbits contribute to conduction. An 
asymmetric potential is therefore needed to observe com- 
mensurability oscillations. This might seem to present 
difficulty, as most devices have symmetric patterns, but 
real structures contain strain that induces an asymmetric 
potential through the piezoelectric effect lia]. The be- 
havior in two dimensions can be explained from the drift 
of the guiding center of cyclotron motion along contours 
of an effective potential, as in one dimension, but the cou- 
pling of motion along x and y causes perturbation theory 
to fail. This coupling also reduces the stability of com- 
mensurability oscillations, and motion becomes chaotic 
at low magnetic fields. Higher Fourier components in 
the potential lead to characteristic signatures in both pxx 
and Pyy near minima of the fundamental oscillations, and 
their detection would verify the theory presented here. 
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